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1. Introduction 

Researchers have increasingly been studying models from economics and from 
the natural sciences where the underlying randomness contains jumps. To give 
an example from financial mathematics, the classical model for a stock price is 
that of a geometric Brownian motion. However, wars, decisions of the Federal 
Reserve and other central banks, and other news can cause the stock price to 
make a sudden shift. To model this, one would like to represent the stock price 
by a process that has jumps. 

This paper is a survey of some aspects of stochastic differential equations 
(SDEs) with jumps. As will quickly become apparent, the theory of SDEs with 
jumps is nowhere near as well developed as the theory of continuous SDEs, and 
in fact, is still in a relatively primitive stage. In my opinion, the field is both 
fertile and important. To encourage readers to undertake research in this area, 
I have mentioned some open problems. 

Section[21is a description of stochastic integration when there are jumps. 
Section |21 describes a number of the most important types of SDEs that involve 
jumps. In Section 0] I discuss a few pathwise uniqueness results, and in Sec- 
tional discuss some results relating to martingale problems. Section|H|concerns 
the regularity of harmonic functions with respect to the corresponding integral 
operators. 

A forthcoming book by Applebaum |Ap| contains much information on 
the subject of this paper, as well as applications to finance. 

"This is an original survey paper. 

t Research partially supported by NSF grant DMS-0244737. 
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2. Stochastic calculus 

Before we can describe SDEs with respect to jump processes, we need to talk a 
bit about the differences between the stochastic calculus for continuous processes 
and for processes with jumps. Some good references for this are the volumes by 
Dellacherie and Meyer DM1 , DM2 , Meyer's course Me , and the books by 
He, Wang, and Yan |HWY| . Jacod [Jd|, Protter [F^, and von Weizsacker and 
Winkler |WW| . I have some notes on the web at |Ba5| : my notes |Ba6] provide 
further background on predictability. 

Suppose we are given a probability space (fi, F, P) and a filtration {.Ft}. 
We suppose that the filtration satisfies the "usual conditions," which means 
that Ft is right continuous and each Ft is complete. We say that X t is right 
continuous with left limits if there exists a null set N so that if u> ^ N, then 
\im u n X u (lu) — Xf(u>) for all t and lim s ^ t X s (ui) exists for all t. The French 
abbreviation for this is "cadlag." Given a process X t that is right continuous 
with left limits, let X t - — lim s | t X s and AX t — X t — X t -. 

A stochastic process X can be viewed as a map from [0, oo) x VL to 
(usually) M. We define a cr-field V on [0, oo) x f2, called the predictable cr-field 
or previsible cr-field, by letting V be the cr-field on [0, oo) x D, generated by the 
class of stochastic processes that are adapted to the filtration {Ft} and have 
left continuous paths. One can show that V is also generated by the class of 
processes H(s,uS) — G(u>)lr a M(s), where G is bounded and F a -measurable. 

Suppose A t is an Ft-adapted process that is right continuous with left 
limits whose paths are increasing, and let us suppose for simplicity that A t is 
integrable for each t. Since A t > A s , then E [A t | F s ] > E [A s | F s ] = A s , so 
trivially A t is a submartingale. By the Doob-Meyer decomposition, there exists 
a predictable increasing adapted process A t such that 

A t = martingale + A t . 

We call A t the compensator or dual predictable projection of A t . If A t has paths 
of bounded variation, we write A t = Af — A^ , where Af and A~[ are increasing 

processes, and then define A t — Af — A^ . 

Before defining the stochastic integral, we discuss the decomposition of 
square integrable martingales. A martingale M t is square integrable if sup t E M t 2 < 
oo. 

If M t is a square integrable martingale that is right continuous with left 
limits, for each integer i let T a = in{{t : \ AM t \ € [2\ 2 4+1 )}, T l2 = inf{t > T a : 
|AM t | S [2 l ,2 I+1 )}, and so on. Since M t is right continuous with left limits, 
Tij — > oo as j ' ^ oo for each i. Let {S, n } be some ordering of the {Tij}. We 
then have a countable sequence of stopping times Si such that every jump of 
M occurs at one of the Si and AMg, is a bounded random variable for each i. 

Let us set = AM Sj l( t > Sj ) and set Mj(i) = Aj( t) - Ai(t). One can 
then prove the following theorem; see |Mej . T.II.ll or Ba5 , Th. 2.3. 

Theorem 2.1 Suppose M t is a square integrable martingale that is right con- 
tinuous with left limits and Mi(t) is defined as above. Then each Mi{t) is a 
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square integrable martingale and Y]°Z 1 Mj(t) converges in L 2 for each t. If 
M t c = M t — J2ili Mi(t), it is possible to find a version of M£ that is a square 
integrable martingale with continuous paths. Moreover, M c and the Mi are mu- 
tually orthogonal. 

The bit about finding a version of M c is due to the fact that the infinite sum 
converges in L 2 , but there is a null set that depends on t. In fact, it is possible 
to arrange matters so that there a single null set. Saying two martingales Ni(t) 
and N2(t) are orthogonal means here that E [Ni(T)N2(S)] = for every pair of 
stopping times S and T. 

Recall that if M t is a martingale, then (M) t is the unique increasing 
predictable process such that M 2 — (M) t is a martingale. The existence and 
uniqueness of (M) t is guaranteed by the Doob-Meyer decomposition. If M t is a 
square integrable martingale whose paths are left continuous with right limits, 
define 

lM] t = {M c } t + J2\&M s \ 2 . 

s<t 

Here M c is the continuous part of M given in Theorem 2.1. One can show, using 
Theorem 2.1, that M 2 — [M] t is a martingale, and in particular, E M 2 = E [M] t . 

With this as background, we can now proceed to a definition of stochastic 
integrals with respect to a square integrable martingale. We want our integrands 
to be predictable. Let us take a moment to explain why this is very natural. Let 
M t — 1 + Pt/\i — t A 1, i.e., a Poisson process minus its mean stopped at time 1, 
and then with 1 added so that M t is nonnegative. Let us suppose that M t is the 
price of a stock, and H s is the number of shares we hold at time s. With this 
investment strategy, it is not hard to see that the net profit (or loss) at time 1 
is Jg 1 H s dM s . Let T be the first time the Poisson process jumps one. If we were 
allowed to choose H s to be zero for s < T and 1 for T < s, our profit at time 
1 would be 1 on the event that (T < 1) and on the event (T > 1); we would 
have made a profit without any risk (if M t has paths of bounded variation, the 
stochastic integral and Lebesgue-Stieltjes integral will coincide). The problem is 
that we looked an instantaneous amount into the future to see when the Poisson 
process jumped. We can't allow that, and the way to prevent this is to require 
H s to be predictable. 

Let us give the definition of stochastic integral. Suppose M t is a square 
integrable martingale with paths that are right continuous with left limits. If 

m 

H s (u) = Y,G i (w)l iaiM] (s), (2.1) 

i=l 

where Gi is bounded and T ai -measurable, define 

N t = H s dM s = G l (M tAbt - M tAOi ). 
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Just as in Ba3 , pp. 43-44, the left hand side will be a martingale and very 
similarly to the proof of |Ba3| . Section 1.5 with [ ] instead of ( ), — [N] t will 
be a martingale, where [N] t — J H 2 d[M] s . 

If if is only P-measurable and satisfies E H 2 d[M] s < oo, approxi- 
mate if by integrands if™ where each if™ is of the form given in i|2.1[l . De- 
fine TV™ = J Q H™dM s . By the same proof as in |Ba3j . Section 1.5, the martin- 
gales N* converge in L 2 . We call the limit N t = f Q H s dM s . The stochastic 
integral is a square integrable martingale, its paths are right continuous with 
left limits, its definition is independent of which sequence if™ we choose, and 
[N] t = f*H*d[M] s . 

We want to generalize the definition of stochastic integral to more general 
processes. For example, a Cauchy process, even if stopped at a fixed time, is 
not square integrable, but we should be able to handle stochastic integrals with 
respect to a Cauchy process by looking at the large jumps separately. 

We say M t is a local martingale if there exist stopping times T n increasing 
to infinity such that for each n the process M t /\T n is a uniformly integrable 
martingale. A semimartingale is a process of the form X t = Xq + M t + A t , where 
Xq is a finite random variable that is .Fo-measurable, M t is a local martingale, 
and A t is a process whose paths have bounded variation on [0, t] for each t. A key 
result concerning semimartingales is the following reduction theorem; see |Me| . 
T.IV.8 or Ba5 , Th. 5.4. This is easy in the case of continuous semimartingales, 
but not at all in the case of semimartingales that have jumps. 

Theorem 2.2 Suppose X t is a semimartingale. There exist stopping times S n f 
oo such that X tA s n = U" + V™, where U n is a square integrable martingale and 
V n is a process whose paths have bounded variation and the total variation of 
V n over the time interval [0, oo] is finite. Moreover, U" = Ug and V t n = Vg 
for t> S n . 

If X t is a local martingale and S n are stopping times such as in Theorem 
2.2, set A7 ASji = {U n Yt for each n and [X] tASn = <* c ) tAS „ + E s <ms„ AX l 
where (U n ) c is defined as in Theorem 2.1. It is not hard to see that these 
definition are independent of the choice of stopping times S n . 

Next we weaken the assumptions on if. A predictable process ff s is 
locally bounded if there exist stopping times R n f oo and constants K n such 
that the process if is bounded by K n on [0, R n }. If if is locally bounded and X 
is a semimartingale, we define J * H s dX s by setting 

H s dX s = I H a ^a n dU^' A g n + / H sA R n dV s n ASn 
Jo Jo 

when t < R n A S n ; the first integral is a stochastic integral with respect to the 
square integrable martingale U n and the second is a Lebesgue-Stieltjes integral. 
Since R n A S n f oo, this defines J H s dX s for all t. It can be shown that this 
integral does not depend on the choice of R n and S n . 

Ito's formula (see [Hi, T.III.3 or [Ba5| . Th. 4.1) is 
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Theorem 2.3 Suppose X is a semimartingale and f is a C 2 function. Then 
f{X t ) is also a semimartingale and we have 

f(X t ) = f(X ) + [ f'(X^)dX s + I / f"(X s _)d(X c ) s 
Jo Jo 

+ J2lf( X °) - /(*-) - f'(Xs-)AX s }. 

s<t 

Note that /'(X s _) is a left continuous process, hence predictable. For ci-dimen- 
sional processes, each of whose components is a semimartingale, the obvious 
generalization holds. 

3. Jump processes 

SDEs with respect to Levy processes. The classical SDE with respect to Brow- 
nian motion is 

dX t = a(X t )dW t + b(X t )dt, (3.1) 

where Wt is a Brownian motion. The simplest analogue of this in the jump case 
is 

dX t = a{X t _)dZ u (3.2) 

where Z t is a Levy process. To be even more specific, we might take Z t to be 
a symmetric stable process of index a. Already, even in this very special case, 
there are several interesting things one can say. Note that we write X t - instead 
of X t in order that the integrand be predictable. 

It would be natural to allow a process to have both a continuous compo- 
nent and a jump component, so one might want to consider the SDE 

dX t = a(X t )dW t + b(X t )dt + a{X t J)dZ t . (3.3) 

Poisson point processes. For many applications, l|3.3[) is not a very good model. 
For example, suppose one wants to model a stock price in such a way that the 
underlying randomness is given by a jump process. For simplicity let us consider 
<|3.2[l instead of l|3.3|l . If Z t has a jump of size z, then X t will have a jump of 
size a(X t -)z. However, one might very well want X t to have a jump whose 
size depends on X t ~ and z, but is not necessarily linear in z. If the underlying 
randomness has a big jump, the behavior of X might be qualitatively different 
from the case where the underlying randomness has a small jump. 

To obtain models with this extra versatility, we need to consider Poisson 
point processes. See |Sk| for details. Let (S, A) be an arbitrary measure space 
(letting S = M. and A be Lebesgue measure will usually do). For each u> e ft 
let /x(w, dt, dz) be a measure on [0, oo) x S. The random measure fj, is a Poisson 
point process if (i) for each set A C S with A(5) < oo the process /x([0, t] x A) is 
a Poisson process with parameter A(^4) and (ii) if A\, . . . , A n are disjoint subsets 
of S with X(Ai) < oo, then the processes /u([0, t] x Aj) are independent. 

Define a non-random measure v by v{[0 : t] X A) = t\(A). If X(A) < oo, 
then /i([0, t] X A) — u([0, t] x A) is the same as a Poisson process minus its mean, 
hence is a martingale. 
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We can define a stochastic integral with respect to the compensated point 
process fj, — v as follows. Suppose H (s, z) — H(uj 1 s, z) is of the form 

n 

H(s, z)=J2 ^M^mMU (*)> ( 3 - 4 ) 
i=i 

where for each i the random variable Ki is J- ai -measurable and A; C S with 
A(ylj) < co. We define 

ft p n 
iV t =/ /ir(s > «)(^-i/)(da,d«)=y)jir < (^-i/)(((a i ,6i]n[0,t])xil < ). 

By linearity it is easy to see that N t is a martingale. It is also easy to see that 
N c = and f 

[N} t = J [ H(s,z) 2 n(ds,dz). (3.5) 
With a little work one can show 



(i 



(N} t = J J H(s,z) 2 v(ds,dz). (3.6) 

Suppose H(s, z) is a predictable process in the following sense: H is mea- 
surable with respect to the a- field generated by all processes of the form l|3.4|) . 
Suppose also that E J Q J H(s, z) 2 v(ds, dz) < oo. If we take processes H n 
of the form l|3.4|l converging to H in an appropriate way, the corresponding 
iV" = J J H n (s,z)d(fL — v) will converge in L 2 , and we call the limit Nt the 
stochastic integral of H with respect to fj,— v. One can show that (|3.5|1 and (|3.6|l 
are still valid. One may think of the stochastic integral as follows: if /i assigns 
mass one to the point (t, z), then N t jumps at this time t and the size of the 
jump is H(t, z). 

Now consider a stochastic differential equation with respect to a com- 
pensated Poisson point process. Look at 

dX t = a{X t )dW t +b{X t )dt + J F(X t -,z)d(vi-v), X = x . (3.7) 

This means 

X t = x + f a(X s )dW s + f b(X s )ds+ f I ' F(X a -,z)(n-v)(ds,dz), 



o 



where Wt is a standard Brownian motion on K. This is a quite general formu- 
lation, as is shown in Cmlar and Jacod |(JJ| . 

Martingale problems. For simplicity we consider here the SDE (|3.7|l with a = 
6 = 0. Suppose / e C 2 and suppose that X t is the solution to 

dX t = J F(X t -,z)d(ji-v). (3.8) 
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By Ito's formula (Theorem 2.3), 

f(X t ) = f(X ) + f f'(X s _)dX s + V[/(X S ) - f(X s _) - f'(X s _)AX s }. 

J " 8<t 

The stochastic integral term is a martingale. We can write f(X s ) as f(X s - + 
AX S ). The jump of X t at time s is equal to F(X S _, z) if n assigns mass one to 
the point (s, z), and so f{X t ) — f(Xo) is equal to a martingale plus 

r [f(X s _ + F(X._,z)) - f(X a _) - f'(X.-)F(X.-,z)]n(d8,dz). 



This in turn is equal to a martingale plus 

[f(X s _ + F(X,_,z)) - /(*,_) - f'{X s _)F{X^,z)}v{dz,ds) 

J [f{X s _ + F(X s _,z)) - f(X s _) - f'(X s _)F(X s _, z)]\{dz)ds. 
If we now set 

Cf(x) = [ [f(x + F(x, z)) - f{x) - f'(x)F(x, z)]X(dz), (3.9) 



we then see that 



f(X t )-f(X )- [ Cf(X a _)ds 
Jo 



is a martingale. Since X has only countably many jumps, then the Lebesgue 
measure of the set of times where X s ^ X s - is 0, and hence f(X t ) — f(Xo) — 
f* Cf(X s )ds = f(X t ) - f(X ) ~ f* Cf(X s ^)ds. Therefore 

f(Xt) - f(X ) - f Cf(X s )ds 
Jo 

is a martingale. 

These calculations (and analogous ones when a or b in (j3.7|l is nonzero) 
are the motivation for what is known as the martingale problem. In this formu- 
lation we take X t to be the canonical coordinate process. One defines an integral 
operator by a formula such as (|3.9() and then says that a probability measure 
V solves the martingale problem started at a point xo if P(Xo = xq) = 1 and 
f{X t ) — f(X ) — J Q Cf(X s )ds is a martingale under P whenever / e C 2 . Note 
here that a solution is a probability, not a process. 

Since the number of large jumps is finite in number over any finite time 
interval and the large jumps do not affect the existence or uniqueness of solu- 
tions, and one wants not to worry about the integrability of J \z\F(x, z)X(dz), 
one often sees operators of the form 

= / [f(x + z)- f(x) - l ([zl < 1) f'(x)z]n(x, dz). (3.10) 
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If n(x, dz) does not depend on x, the reader will note that in this case (|3.1U|) is 
the infinitesimal generator of a Levy process. 

An example of an operator such as l|3.10fl is to let 

n{x,dz) = ^ l+a dz, (3. II) 

where a is bounded above and below by positive constants, and < a < 2. If 
a is constant, we have the infinitesimal generator of a symmetric stable process 
of index a. 

The extension of the notion of martingale problem to higher dimensions 
causes no problems. One also often sees operators of the form C + Ai, where C 
is as in <|3.l0fl and 

Mf(x) = \a{xff"(x) + b(x)f'(x). 

This corresponds to adding a diffusion and drift term to the jump operator. 
Again, the higher dimensional analogues are what one expects. 

Pseudodifferential operators. With C as in I|3.l0ll . let us see what we get when 
we apply C to the function f(x) — e mx . We obtain 



Cf{x) = e mx J [e mz - l - iuzl { \ z \< x) )n{x, dz). 

In the case where n(x,dz) does not depend on x, i.e., n(x,dz) — n(dz) for all 
x, this implies that the Fourier transform of C is ip(u), where 



ip(u) = J [e luz - 1 - mzl(| z |< 1) ]n(dz). 

More generally, when n{x, dz) does depend on x, we let 

ip(x, u) — J [e" tz — 1 — mzl(| 2 |< 1 )]n(a;, dz), 

and we call ip{x,u) the symbol corresponding to the operator C. 

By the uniqueness of the Fourier transform, we can specify an operator by 
presenting n(x, dz) or we can specify the operator by giving its symbol ip(x, u). 

Dirichlet forms. Suppose one considers the operator Cf(x) — (a(x) f (x))' . For 
this to make sense in terms of the usual notion of derivative, we need a to be 
differentiable. If g is C°° with compact support, an integration by parts gives 

g(x)Cf(x)dx = / g'(x)a(x)f'(x)dx. 



The expression on the right makes sense for any measurable and bounded a, 
as long as / and g are differentiable with compact support, say, and is called a 
Dirichlet form, written E{f,g). 
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Under minimal hypotheses, a Dirichlct form will determine the process. 
See |FOT| for the basic theory of Dirichlet forms. 
If we define 

£(f, f) = J J lf(y) - f(x)] 2 J(dx, dy), (3.12) 

where the measure J is symmetric, the associated process will be of jump type. 
A special case is when J(dxdy) = \x— y\~ 1 ~ a dx dy, and the associated process is 
a symmetric stable processes of index a. We will briefly consider some properties 
of processes associated to the Dirichlet form given in (|3.12f) . called symmetric 
jump processes, in the case where 

J(dx,dy)= a [ X ^ +a dxdy (3.13) 

and also higher dimensional analogues, where < a < 2, a is bounded above 
and below by positive constants, and a[x, y) — a(y, x) for all x and y. Here the 
questions concern not uniqueness, since that is covered by the general theory of 
Dirichlct forms, but instead properties of the corresponding process. (But see 
also US]-) 

4. Pathwise uniqueness 

When is the solution to (|3.2|l or l|3.7|l pathwise unique? When the coefficients 
are Lipschitz coefficients, then the standard Picard iteration procedure proves 
uniqueness. For example, for (|3.7|l we have the following, due to |Skj . 

Theorem 4.1 If a and b are bounded and Lipschitz, Jsup x \F(x, z)\ 2 \(dz) < 
oo and 

' \F(x,z)-F(y,z)\ 2 \(dz)<c 1 \x-y\ 2 



for all x, y, then there exists a solution to \3. r fy and that solution is pathwise 
unique. 

Concerning 1)3.2(1 . again, if a is bounded and Lipschitz continuous, one 
would expect pathwise uniqueness, and indeed that is the case. In one dimension, 
however, in view of the result of Yamada-Watanabe [YWj for SDEs with respect 
to a one dimensional Brownian motion, one would hope that much weaker condi- 
tions would suffice for uniqueness. For example, the Yamada-Watanabe sufficient 
condition for pathwise uniqueness for diffusions is that Jg p(x)~ 2 dx = oo for all 
e > 0, where p is the modulus of continuity: \o~(x) — o~(y)\ < p(\x — y\) for all x 
and y. 

For solutions to l|3.2|l we have the following; see |Km2j and |Ba4| . 

Theorem 4.2 Suppose a > 1, Z t is a one dimensional symmetric stable process 
of index a, and \a(x) — a(y)\ < p(\x — y\) for all x and y and 

, 7 dx = oo 
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for every e > 0. Then for every xq G R, the solution to 
dX t = a(X t -)dZ t , X = x , 

is pathwise unique. 

This result is sharp. As a corollary, if a is Holder continuous of order 1/a, 
then we have pathwise uniqueness. Note that the smaller a, the more continuity 
is required. 

What happens when a < 1? When a < 1, the paths of Z t are of bounded 
variation on finite time intervals, and I initially thought only continuity and 
appropriate boundedness of a sufficed |Ba4j . However the solution constructed 
there, although measurable with respect to the c-field generated by Z, is not 
adapted. In fact, when a < 1, the condition that a be Lipschitz continuous 
turns out to be sharp; see |BBC| . Probably the Lipschitz continuity condition 
can be weakened by a logarithm term, but it does not suffice to let a be Holder 
continuous of order (3 for any (3 < 1 . 

For diffusions Nakao |Na) showed that bounded variation of a suffices 
for %\.\\ to have a unique solution, and LeGall [LO| improved this to a hav- 
ing finite quadratic variation. What is the appropriate analogue of this for the 
equation (|3.2(l ? 

A paper concerning pathwise existence along somewhat different lines is 
Williams jWij . Here the stochastic integral is not the Ito integral, but instead 
the rough integral of Lyons (cf. |Ly| ) , and the issue is existence. 

Other papers related to pathwise properties of SDEs with jumps are |Roj 
on backward SDEs, |JMW) on approximations of solutions, and |AT| . |Fu| . | FK| . 
and |Kuj on flows of solutions. 

5. Martingale problems 

In this section let us discuss existence and solutions to martingale problems such 
as those mentioned in Section [3J When there exists a non-degenerate diffusion 
component, that term dominates the situation, and uniqueness holds. We have 
the following theorem of Stroock St 

Theorem 5.1 Suppose 

Cf(x) = a l3 {x)^l-(x) +j^h l {x)^-(x) (5.1) 

+ 1 [f(x + z)-f(x)-l Qzl < 1) Vf(x)-z]n(x,dz) (5.2) 

JR d \{0} 

and 

(a) a is bounded and continuous and strictly elliptic; 

(b) b is bounded and measurable; 

(c) j A i+\ z \2 n(x, dz) is bounded and continuous for each A C M d \ {0}. 
Then there is a unique solution to the martingale problem for C started at any 
xq e R d . 
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This result should also be compared with that of Komatsu |Kml| and that of 
Lepeltier and Marchal |LM| . 

Suppose now that there is no diffusion component present. For early 
works see |Tsl| . |Ts2j . and |TTW] . A special case of a result of Komatsu |Km4j 
is the following. 

Theorem 5.2 Suppose 
Cf{x) = [ [f(x + z) - f{x) - l (k |<DV/(x) • z] 

jR d \{0} 

for some a S (0,2), n*(dz) = swp x n(x,dz), and 

J (1 A \z\ a )n* (dz) < oo. 

Then there is a unique solution to the martingale problem for C started at any 
point xq. 

See also [NT] . 

Martingale problems are closely related to the notion of weak uniqueness 
for SDEs. Concerning weak existence of solutions for (|3.2|1 see jPZj . |Zalj . and 
|7a2] . 

Another approach that has been explored has been the use of pseudod- 
iffcrential operators. See works of Komatsu, Jacob, Hoh, and Schilling: Km3 , 

EH!, EEJ, [EEL]. GE2I, GEH], GEU, CP, EEU, PS3, EE3, GE1, Call, and 

[Ja5] . 

A recent result along these lines is that of Kolokoltsov K12 . If 
i = are the symbols associated to Levy processes satisfying certain 

conditions and a^(x) are suitable functions of x, then there is a unique solution 
to the martingale problem associated to the symbol 

n 

ipix, u) = aiix)ipiiu). 
t=i 

Can one obtain a similar result when we have an infinite sum instead of a finite 
one? 

In most of the above results, the operators are the sum of a finite number 
of integral terms plus a lower order term. There are fewer papers that handle 
variable order terms without assuming a considerable amount of continuity in 
the x variable. One such paper is |Bal| . This paper deals with one-dimensional 
processes, and the exact conditions are rather complicated. One can give an 
example, though, that illustrates the theorem. Suppose C is given by (|3.10() . 
where 

n(x,dz) = — , . dz. (5.3) 

v ' ' \ z \l+a(x) y ' 

Qualitatively, this says that at a point x, the process behaves like a symmetric 
stable process of order aix). Note that the order depends on x. For this example, 
the conditions of |Bal reduce to 



\d+a 



dz + noix, dz) 
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Theorem 5.3 Suppose a{x) is Dini continuous, bounded above by a constant 
less than 2 and bounded below by a constant greater than 0. If C is given by ^3.1(J\) 
with n{x, dz) defined by then there exists a unique solution to the martin- 

gale problem for C started at any point io el. 

Dini continuity of a means that there exists p such that \a(x) — oc(y)\ < 
p(\x — y\) for all x and y and 

^ - dx < oo 



o 

for all e > 0. 

Other papers that deal with properties of processes associated to oper- 
ators of variable order without assuming a great deal of smoothness include 

H52I, lEi, HiU, IKS3, Em, El 

It would be greatly desirable to have a uniqueness theorem for variable 
order operators where the hypotheses are relatively simple. 

6. Regularity of harmonic functions 

In this section we give some regularity results for harmonic functions related to 
the integral operators C such as those similar to 1^.111) . Before doing this, we 
mention that there has been a great deal of recent work on the potential theory 
of symmetric stable processes in domains, such as Fatou theorems, boundary 
Harnack principle, intrinsic ultracontractivity, Green function estimates, etc. See 
jChj for a survey of some of the recent work. Little of this has been extended to 
more general jump processes, but it would be worthwhile doing so. 

First we consider the Harnack inequality. Suppose Ch{x) = in a domain 
D C K d . Here the operators C that we are considering are of the form 



»(*>*) = (6-2) 



Cf(x)= [f(x + z)-f(x)-l {lzl < 1) Vf(x)-z]n(x,z)dz, (6.1) 

jR d \{0} 

where 

a(x, z) 

w 

a G (0,2), and a is bounded above and below by positive constants. It may 
so happen that what one would like to consider as a harmonic function is not 
regular enough to be in the domain of C. So a more general definition of harmonic 
function is to say that h is harmonic if h(X tATD ) is a martingale. Here td — 
inf{i : X t D}, the first exit time of the domain D, and X t is the Markov 
process associated to the operator L. We have the following Harnack inequality 

EH]. 

Theorem 6.1 Suppose C is given by ]6.1)) . n by n(x,z) — n(x,—z) for 

all x and z, xq £ M. d , R > 0, h is nonnegative and harmonic in a ball B(xq, 2R) 
of radius 2R centered at xq, and h is bounded in R d . Then there exists a constant 
C\ not depending on h, Xq, or R such that 

h{x) < cih(y), x, y € B(x Q ,R). 
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This theorem has been extended by |SV| to more general n, but still ones 
that are essentially stable of some order a. 

One of the reasons one is interested in Harnack inequalities is that they 
often imply that harmonic functions must be continuous. In BL1 it is also 
proved that 

Theorem 6.2 Under the same conditions as Theorem 6.1 there exist constants 
C2 and 7 £ (0, 1) such that 

Hx)~h{y)\<c 2 {^ 1 ^Y\\h\\ 00 . 

In other words, harmonic functions are Holder continuous. 

What if the function n is not of the form (|6.2(l ? Suppose there exist 
constants C3 , C4 and < a < (3 < 2 such that 

~ n{x ' v) - \x-l\^ l*-^ 1 ' ( 6 - 3 ) 

and with some appropriate conditions on n(x, y) when \x — y\ > 1? In |BKj it 
was proved that the conclusion of Theorem 6.1 holds provided that (3 — a < 1 
and one allows c\ to depend on R. It was shown that the dependence on R 
cannot be dispensed with. 

It is still unknown whether a Harnack inequality must hold in the case 
/3 — a > 1, or whether a counterexample exists. 

If one turns to symmetric jump processes, that is, ones determined by 
the Dirichlet form H3.12[l . that the analogue of Theorem 6.1 holds is implicitly 
proved in |BL2| ; see also |CK| . In the paper |BL2| upper and lower bounds on 
the transition probability densities were found, extending and improving results 
of Komatsu |Km5| . 

I believe there is also an analogue of jBKj . but that in this case the full 
range 0<a</3<2is allowed. Whether the analogue of Theorem 6.2 holds in 
this full range is at present unknown. 

If one assumes more regularity in the operator £, that is, stronger con- 
ditions on n, then more can be said about regularity. See |K12| . 
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